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(N ! Abstract 

This work presents a method for the ab initio calculation of isotope shift in atoms and ions with 
one valence electron above closed shells. As a zero approximation we use relativistic Hartree-Fock 
and then calculate correlation corrections. The main motivation for developing the method comes 
, from the need to analyse whether different isotope abundances in early universe can contribute 

to the observed anomalies in quasar absorption spectra. The current best explanation for these 
anomalies is the assumption that the fine structure constant a was smaller at early epoch. We test 
the isotope shift method by comparing the calculated and experimental isotope shift for the alkali 
and alkali-like atoms Na, Mgll, K, Call and Ball. The agreement is found to be good. We then 
calculate the isotope shift for some astronomically relevant transitions in Sill and SilV, Mgll, Znll 
and Gell. 

*Q I PACS numbers: 31.30.Gs, 31.15.Md, 31.25.Jf 
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I. INTRODUCTION 



Recent studies of quasar absorption spectra reveal a possible change in a since the early 
universe One of the possible major sources of systematic effects in these studies is that 
the isotopic abundance ratios in gas clouds in the early universe could be different to those 
on Earth. This may provide an alternative explanation for the observed variation in spectra 
0. In order to test this possibility it is necessary to have accurate values for the isotope 
shift (IS) for the relevant atomic transitions. Experimental data is available for only a very 
few of them. Therefore, accurate calculations are needed to make the most comprehensive 
analysis possible. 

Another motivation for accurate isotope shift calculations comes from the possibility to 
use isotope shift to study atomic nuclei. Comparing calculated and measured IS allows one 
to find the change in nuclear charge distribution from one isotope to another. Studying IS for 
these purposes has long history, for a good review and tables see e.g. or for more recent 
values in many atoms and ions j^|. However any progress in the accuracy of calculations is 
of significant importance. 

Isotope shift is also important in searching for small, charged black holes. The mass deficit 
of the observable universe (dark matter) may be explained by supposing the existence of 
non-vanishing "elementary" black holes of the Planck mass. Such black holes may have 
an electric charge, giving rise to the possibility of an atom made from electrons orbiting 
a positively charged black hole instead of a nucleus. Such black hole atoms would have 
spectra shifted with respect to usual nuclear atoms due to the effectively infinite mass and 
zero volume of the nucleus p • This shift is simply an extreme example of the regular isotope 
shift, and can be calculated using the method described in this paper. A search for these 
spectra would verify the existence of elementary black holes, or any other stable, very heavy 
particles, e.g. "strange matter" . 

In this paper, we develop a method for calculating the isotopic shift of atoms and ions 
that can be treated as a single electron above a closed shell core. These include the alkali 
metals, as well as other atoms and ions such as Znll that are fairly well approximated by a 
single electron above a closed subshell core. 

We test our technique by calculating IS for light alkaline atoms as well as for a relatively 
heavy ion, Ball. Since different contributions dominate in light and heavy atoms and good 
agreement with experiment has been obtained in both cases, one can confidently say that 
the technique has been reliably tested. We then apply the technique to calculate isotope 
shift for astronomically relevant transitions in Mgll, Znll, Gell, Sill and SilV. 

II. METHOD 

The isotope shifts of atomic transition frequencies come from two sources: the fi- 
nite size of the nuclear charge distribution (the "volume" or "field" shift), and the fi- 
nite mass of the nucleus (see, e.g. 0). The energy shift due to recoil of the nucleus is 
(l/2M)p^r = (l/2M)(Spj) 2 . Furthermore this "mass shift" is traditionally divided into 
the normal mass shift and the specific mass shift (SMS). The normal mass shift is given by 
the operator (l/2M)T,pf, which is easily calculated from the transition frequency. The SMS 
operator is (l/M)S i<: ,pj ■ pj which is difficult to evaluate accurately. 

The shift in energy of any transition in an isotope with mass number A' with respect to 
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an isotope with mass number A can be expressed as 



5v 



A' A 



(k NMS + ksMs) ( j t - \) + F5(rY> A , (1) 



where the normal mass shift constant is 

kNMS = ~ 182^888 (2) 

and (r 2 ) is the mean square nuclear radius. In this paper we develop a method for calculating 
the specific mass shift and field shift constants, ksMs and F respectively. It is worth noting 
that in this paper we use the convention bv A ' A = v A — v A . 

To determine spectral shifts in black hole atoms relative to normal atoms, one must take 
A 1 — > oo, and because the black hole has effectively zero volume 5(r 2 ) = — (r 2 ) = — |i? 2 uc . 



A. Specific mass shift 

It is well known that calculations using many-body perturbation theory (MBPT) in the 
residual Coulomb interaction give very poor convergence for atoms with many electrons. 
Therefore all order techniques were developed in earlier works to calculate energy levels, 
transition amplitudes, etc. (see, e.g. [E IS HI IHJ)- ^ i s natural to expect that an all order 
technique is needed for the isotope shift as well. In the recent work by Safronova and Johnson 
(Ref. the SMS was calculated to third order of MBPT. Their results verify that there 
is very poor convergence for this operator. A first or second order calculation cannot even 
guarantee that the sign of the final result will be correct. A third order calculation gives the 
correct sign, however it is unknown how much fourth and higher order diagrams contribute. 

Our method is to include the specific mass shift directly into an energy calculation from 
the very beginning. The SMS is a two-body operator prP2 and can be added to the Coulomb 
potential, Q = + Api • P2 (see Appendix,). The operator Q replaces the Coulomb 

operator everywhere that it appears in an "all order" energy calculation that includes certain 
chains of diagrams to all orders. We vary the isotope shift scaling factor A and calculate the 
level energies E. The gradient of an E vs. A graph then gives us the SMS matrix element, 
usually denoted k S MS- 

As a zero approximation we use relativistic Hartree-Fock (Hartree-Fock-Dirac) method. 
At the Hartree-Fock stage we include the isotope shift with the exchange potential and iterate 
to obtain self-consistent "dressed" wavef unctions. The SMS matrix element at this stage 
is roughly equivalent to the first-order and second-order contributions of the one-particle 
operator, plus higher-order corrections from the random-phase approximation in Ref. 
(pC 1 ) 4- 5"( 2 ) in their notation). It is interesting to note that ksMS at this stage does not give 
meaningful results, and can even be of the wrong sign (it is labelled as "HF" in Table |1J). 
We need to include correlation corrections in order to obtain any reasonable accuracy. 

We include correlation effects to second order of MBPT; this approach gives good accu- 
racy for energies. The perturbation is the difference between the exact and Hartree-Fock 
Hamiltonians V = H — -£/hf- The many-body corrections start in second order; all first 
order diagrams have been included in the self-consistent Hartree-Fock procedure. There are 
just four second order diagrams, illustrated for example in Ref. Q. However now we replace 
the Coulomb operator in these diagrams with our two body operator Q, and use dressed 
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TABLE I: Comparison of specific mass shift constants, kgMSi obtained by including various 
contributions in the energy calculation (all in GHz.amu). Our values "HF" are approximately 
equivalent to the values labelled "pW + ,S'( 2 )" here and in Ref. 11 1. 



Ion State 



Ref. [Ill 



This work 
pW + Final HF Final 



Nal 3s 


-54 


54 


-52 


69 


3Pl/2 


-67 


-43 


-68 


-40 


3 P3/2 


-67 


-43 


-67 


-39 


Mgll 3s 




38 


-171 


83 


SPl/2 




-324 


-408 


-296 


3 P3/2 




-323 


-402 


-290 



Hartree-Fock wavefunctions as described above. This gives our final value of ksus as listed 
in Table |U 

An estimate of the size of neglected higher order diagrams can be obtained by calculating 
a new value of ksMS that includes some higher order chains of diagrams. We create an 
operator E for the second order correlation effects in each partial wave, defined by 



5E^ = (n\E\n) 



(3) 



This "sigma" operator is then added to the exchange potential in the Hartree-Fock Hamilto- 
nian for the valence electron, H = /7hf + E. Thus we calculate the single electron Brueckner 
orbital. We also include a scaling factor, /, with this sigma operator (E — > /E) in order to 
fit the experimental energy. The rescaling of E simulates some higher order correlation cor- 
rections that were omitted in our calculation. Including E in our Hartree-Fock calculation 
allows us to obtain another value for ksMS, an d the difference between this new value and 
the old second-order value gives us an estimate of the error. 

In Table H] we compare our results to those obtained by Safronova and Johnson (Ref. 
who calculated the isotope shift to third order in standard MBPT. Our values are signifi- 
cantly different, in particular the Mgll 3s level shift is more than twice that given in Ref. 
and also which are broadly in agreement with each other. However we find that this 
difference is not so large in the 3p — 3s transition due to cancellations of higher order terms 



between the two levels. We compare our final results and those of Ref. |ll[ with experiment 



in Table HT1 Agreement is at the level of 1%, which is much smaller than our error estimates. 

Our method includes many-body diagrams that are different to Ref. 11], including some 
chains of diagrams in all orders. It is also arguably simpler to implement since it doesn't 
require term-by-term calculation of a large number of matrix elements. 



B. Field shift 

Our method for calculating the field shift (also called the volume shift) is similar to that 
used for the SMS calculation. We add a perturbation due to the field shift to the nuclear 
potential, and then calculate the energy directly. The field shift potential is 

5U(r) = A ( U(R + 6R, r) - U(R, r) ) (4) 
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TABLE II: Comparison of the specific mass shift of transitions in Na and Mgll with experiment. 



Isotopes Transition SMS (MHz) 

Ref. [11] This work Exp. a 
23 " 22 Na 3p 1/2 - 3s 192 214(48) 215(1) 

214(2) 

3p 3/2 - 3s 192 212(48) 214 
26-24 MgII 3p3/2 _ 3s n57 1196 ( 18 ) 1207(100) 

"these values are extracted from IS experiments in Table IvTl 



where R is the nuclear radius and U(R, r) is the nuclear potential. To obtain the change 
in energy of a state due to the field shift, we add this potential to the nuclear potential in 
our Hartree-Fock calculations. By introducing a scaling factor A we can check linearity and 
increase the size of the effect. To calculate the field shift constant, we take the gradient of 
a level energy E vs. A graph and extract F as 

F= SEy_ 5 1 dE 

S(r 2 ) 32R5RdX ' U 

Note that the field shift constant, F, is defined here as having opposite sign to the definition 
in some previous papers, e.g. Ref. and [IsT ]. This equation assumes that the nuclear 
charge distribution can be approximated as a uniformly charged sphere of radius R. 

Higher order correlation effects are smaller in the field shift than in the specific mass 
shift, and are localised at the nucleus. We can include them easily by creating Brueckner 
orbitals (solutions of the modified Hamiltonian, H = Hup + E) using a second order sigma 
operator, defined in Equation El The field shift is proportional in first order to the square 
of the wavefunction at the nucleus. Hence we can include higher order effects quite simply 
by multiplying the matrix element by the square of the ratio of the Brueckner wavefunction 
to the Hartree-Fock wavefunction at the nucleus, 

Rrf+t, 

We tested this in Ball, as well as in lighter atoms, because in barium the field shift 
dominates over the mass shift, and there is a lot of experimental data to compare with. 
Our results for field shift constants in Ball were found to be consistent within a few percent 
of the previous theoretical work by Martensson-Pendrill (Table ITTT| . In the same paper, 
Ref. [l3], they say that they have underestimated the 6p 1//2 field shift constant by around 7%, 
based on the difference between experimental and theoretical calculations of the hyperfine 
A constant. Also, their 6s constant is said to be overestimated, leading to a corrected value 
of F 6p3/2 _ 6s = —4.20(13) GHz/fm 2 . These corrected values are in better agreement with the 
ab initio values obtained in this work (F 6p3/2 _ 6s = —4.076 GHz/fm 2 ). 

Using King plots 0] we can extract the ratios of field shift constants for different transi- 
tions, provided we have experimental data for a number of different isotopes. In Table IVI we 
compare our calculated values of these ratios with those obtained by combining the data in 
several different experiments and transitions. Our values were found to be consistent with 
this experiment to within 5%. We have used a simple weighted least squares fit to obtain an 



-0; 



HF+S 



(0) 



^hf(O) 



(6) 
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TABLE III: Level field shifts in Ball states. Due to a difference in the definition of F, the values 
calculated in Ref. 0] have been presented here with opposite sign. 



State 


Energy 


F (MHz/fm 2 ) 




(cm -1 ) 


Ref. [13] This 


work 


6s 


-80686.87 


4096 


3851 


6 Pl/2 


-60425.31 


-111.1 


150.1 


6P3/2 


-58734.45 


-242.6 


225.4 


5d 3/2 


-75813.02 




-1223 


5d 5/2 


-75012.05 




-1148 



TABLE IV: Calculated SMS and field shift constants in Ball transitions. 



Transition 


Wavelength (nm) k$MS 


(GHz.amu) F (MHz/fm 2 ) 


6P1/2 - 6s 


493 


105 


-4001 


6P3/2 - 6s 


455 


257 


-4077 


5d 5/2 - 6s 


1762 


-550 


-4999 


6pi/ 2 - 5d 3/2 


650 


653 


1073 


6P3/2 - 5d 3/2 


585 


805 


997 


6P3/2 - 5d 5/2 


614 


807 


922 



experimental value for the ratio F^/F^. A two point formula was used for ratios involving 
the 1762 nm transition as experimental data exists only for the Su 134,138 and Sv 13e ' 138 j so ^ p e 
shifts 0. Other ratios were extracted by the groups that performed the experiments, with 
much higher accuracy. 



TABLE V: Ratios of field shift constants in Ball states. In the second column we list the measured 
values, obtained using King plots. In some cases we obtain the ratio ourselves by combining the 
results of two separate studies. 



Transitions 






FJFp 


( <*/P ) 




This work 


Experiment 


6pi/ 2 - 6s / 6p 3/2 


— 6s 


0.982 


0.975(3) a 


6p 3 /2 - 6s / 6p 3/2 


- 5d 5/2 


-4.42 


-4.50(6) 6c 


6P3/2 - 6s / 5d 5/2 


— 6s 


0.816 


0.82(4) M 


5d 5/2 - 6s / 6p 3/2 


- 5d 5/2 


-5.42 


-5.5(3) cd 


6Pi/ 2 - 5d 3/2 / 6p 3/2 


- 5d 3/2 


1.076 


1.087(4) c 


6P3/2 - 5d 5/2 / 6p 3/2 


- 5d 3/2 


0.925 


0.961(3) c 


6P3/2 - 5d 3 /2 / 6p 3/2 


— 6s 


-0.245 


-0.2312(6) c 



a Wendt et al. 1^ 
b Wendt et al. [ia_ 
c Villemoes et raM17l| 
d Zhao et al. [li 
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TABLE VI: Comparison of experimental values of the specific mass shift with calculated theoretical 
values. The experimental values were extracted by subtracting the NMS and field shift (FS) from 
the experimental IS. 



Isotopes 


Transition 


Energy 
cm 


IS (exp) NMS FS 
(MHz) (MHz) (MHz) 


SMS (MHz) 
This work Exp. 


Z6 — ZZ~\j 

in a 


3P1/2 - 


3s 


loyoo. lo 


( DO.O{ 1 ) 


ooi 


Q 

-o 


214(48) 


215(1) 








( oo.y^i.y ) 








214(2) 




3 P3/2 _ 


3s 


±oy i o.oo 


(O I . ( ZyZQ ) 




Q 
-O 


212(48) 


214 


26 — 24iv/r„TT 
IVLgll 


3P3/2 - 


3s 


35/oO.y / 


3050(100) 


1185 


-42 


1196(18) 1207(100) 


41 — 39 


4 Pl/2 - 


4s 


12985.17 


235.25(75) 


267 


-13(5) 


-32(21) 


-19(6) 




4d l/2 ~ 


■J.S 


27398.11 


K Q K I O ^ f 

ooo^y y 


004 


-13(5) 


20(30) 


34(13) 


A'i — 40 /~\ TT 

^ * u Call 


4 Pi/2 - 


4s 


25191.54 


706(42)9 


723 


-36(3) 


22(1) 


19(45) 








672(9) h 








-15(11) 










685(36)' 








-2(39) 




4 P3/2 - 


4s 


25414.43 


713(31) 9 


729 


-36(3) 


-5(1) 


20(34) 








677(19) fe 








-16(22) 










685(36) s 








-8(39) 




3d 3 / 2 - 


4s 


13650.21 


4180(48)9 


392 


-47(4) 


3502(217) 


3835(52) 




3d 5/2 - 


4s 


13710.90 


4129(10)9 


393 


-47(4) 


3487(215) 


3783(14) 




4 Pl/2 ~ 


3d 3 / 2 


11541.33 


-3464.3(3.0^ 


331 


12(1) 


-3479(218) 


-3807(4) 






-3483(40)9 








-3826(41) 




4 P3/2 - 


3d 3/2 


11764.22 


-3462.4(2.6)^ 


337 


12(1) 


-3507(217) 


-3811(4) 






-3446(20)9 








-3795(21) 




4 P3/2 - 


3d3/2 


11703.53 


-3465. 4(3. 7) J 


336 


12(1) 


-3492(216) 


-3813(5) 






-3427(33)9 








-3774(34) 


1 1 

a Pescht et al. ]u| 
6 Huber et al. 0,— ■ 
c Gangrsky et al. 
d Drullinger et al. XL3 
e Touchard et al. \9.M 
^Horback et n l 
9 Kurth et al. 25] 
''extracted from Marti 
'Maleki and Goble[23 
J N6rtershauser et al. 


onsson- 


Pcndrill et al. 


□ 











C. Alkaline ions 

We compare our results with experimental data for alkaline ions in Table IVII In alka- 
line ions it is more valuable to compare only the specific mass shift with those extracted 
from experiment, than to compare the entire isotope shift. This is because the mass shift 
dominates strongly in these ions, and also because the SMS is generally considered more 
difficult to calculate. We have removed the field shift and the normal mass shift from the 
experimental values of the isotopic shift in order to obtain an experimental value for the 
specific mass shift. The field shift values used in Table EU were calculated using the above 
method. While our calculation of F has been shown to be good, the field shift also depends 
on having knowledge of 5(r 2 ) for the relevant isotopes. 

For Na we use the value quoted in Ref. [TJ of <5(r 2 ) 23 - 22 = 0.20 5(3) fm 2 . This value 
is only from an empirical fit, and shouldn't be trusted too far. The field shift is very 
small in this atom, so the errors don't matter too much. For Mgll we have the used value 
#(r 2 ) 26 ' 24 = 0.55 fm 2 from another empirical fit, the equation -R n uc = 1.1A 1 / 3 fm. This is very 
poor, but in this case the field shift is small even in relation to the error in the experimental 
isotope shift. In Table IVT1 we have not included an error contribution for the field shift in 
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TABLE VII: Mass and field shift constants for some useful transitions. 



Ion 


Transition 


F 


Mb 


ksMS 






(MHz/fm 2 ) (GH 




( CtTTz anil]) 


Mgll 3 Pl/2 - 3s 


-127 


-587 


-373(42") 




3 P3/2 - 3s 


-127 


-588 


-373(6) 


Sill 


4s - 3p 1/2 


171 


-1077 


1257(291 

-LiJL/ 1 1 til/ 1 




4s - 3p 3/2 


171 


-1072 


1243(28) 


SilV 


3pi/ 2 - 3s 


-484 


-1172 


-1535(11) 




3P3/2 ~~ 3s 


-485 


-1180 


-1505(7) 


Znll 


4pi/ 2 - 4s 


-1596 


-797 


-1310(69) 




4p 3 /2 - 4s 


-1596 


-812 


-1266(69) 


Gell 


5s - 4p 1/2 


1088 


-1026 


1046(69) 




5s - 4p 3/2 


1083 


-997 


960(62) 



either of these atoms, since we really don't know how accurate these approximations are. 

The values of 5(r 2 ) are known for K and Call from muonic x-ray experiments, allowing 
us to calculate the field shift much more accurately. This is fortunate because the SMS is 
relatively small for the p — s transitions in these atoms, and hence the field shift plays a 
much larger role. We use the values <5( r 2 ) 41 ' 39 = 0.117(40) fm 2 for K from Ref. @, and 



<5(r 2 ) 43,40 = 0.1254(32) fm 2 for Call from Ref. |30|. In Call the change in mean square 
nuclear radius is given to high precision, so we have included an additional error of 5% in 
the field shift that comes from the constant F. This is a pessimistic estimate of error based 
on the accuracy we achieved calculating F for transitions in Ball. 

Table IVT1 shows that our method can reliably calculate the isotope shift in alkaline atoms, 
including those transitions with a large specific mass shift. 



III. RESULTS 

We have shown that our method works in atoms for which we have available experimental 
data (Section inj- In Table IvTTl we tabulate values for the mass and field shift constants for 
some astronomically useful transitions. We have not given errors for F, however we can say 
that they are less than 5% based on comparison of calculation with experiment in Ball. 

In Table IVTTT1 we present the results of isotope shift calculations between common isotopes 
of astronomically important ions. We have used the IS constants presented in Table IVTT1 with 
Equation ^ in order to calculate the isotope shift between particular isotopes. Just before 
submission of this paper, results of measurements for Znll were brought to our attention 
|3l| . These results matched our prediction extremely well. 

In Gell and Sill, the specific mass shift cancels the normal mass shift entirely, making 
the field shift, and hence S(r 2 ), important. We have just used the empirical formula R nnc — 
1.1A 1 ' 3 fm 2 to obtain values of 5(r 2 ) in these ions. This is extremely rough, and although 
it seems to work in Znll to within a 20% accuracy based on the experimental data given, 
we really don't know if this holds for Si and Ge at all. In the SilV transitions presented, 
it is less important to have good values for the field shift because there is no cancellation 
between the NMS and SMS. 



8 



TABLE VIII: Summary of isotope shift values for astronomically relevant alkali-like ions. The 
experimental value for the 4p 3 / 2 — 4s transition in Znll is 676(6) MHz as quoted in 3l(. We have 
presented two errors, the first is our uncertainty in k$MS an d the second is the uncertainty in our 
field shift, which is mainly due to lack of knowledge of 5{r 2 ). A negative shift means that the sign 
is opposite to the normal mass shift. 



Isotopes Transition Energy Isotope Shift 

(cm- 1 ) (MHz) 
3 °- 28 SiII 4s-3p 1/2 65495.1 -375(70)(11) 

4s-3p 3/2 65208.1 -351(67)(11) 
30-28 Snv 3p 1/2 -3s 71289.6 6294(26) (31) 

3p 3/2 - 3s 71749.9 6241(18)(31) 
66 " 64 ZnII 4p 1/2 - 4s 48480.6 653(32) (78) 

4p 3/2 - 4s 49354.4 632(33) (79) 
74 - 70 GeII 5s-4p 1/2 62402.4 491(53)(101) 
5s-4p 3/2 60635.3 533(48)(101) 



IV. CONCLUSION 

We have presented a method for the calculation of the isotope shift in atoms and ions 
that can be approximated as having one valence electron above a closed shell. Our results 
are shown to be in good agreement with isotope shift experiments in both light and heavy 
atoms, which are good tests for the mass shift and field shift respectively. 

We have used the method to predict values of the isotope shift in astronomically relevant 
transitions in Sill and SilV, Mgll, Znll and Gell. Recent experiments measured the isotope 
shift of a p — s transition in Znll, and the results were in excellent agreement with our 
prediction. These values are needed in order to examine systematic effects in observations 
of spectral line shifts in quasar absorption spectra that suggest a variation in a 0. Our 
calculations could help provide another explanation for the observed shifts, depending on 
whether isotopic abundances were different in gas clouds in the early universe. Alternatively, 
our calculations may strengthen the arguments in support of a varying a. 

Further work needs to be done in order to obtain the isotope shift for more complex 
atoms with more than one electron in their outer shell. The general method of including 
the isotope shift operator with the two-body Coulomb operator may be used in such cases, 
although the energy calculation itself is more complicated. 
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APPENDIX: MATRIX ELEMENT OF THE TWO-BODY OPERATOR 



The two-body operator used is this work is the sum of the Coulomb interaction operator 
and the "rescaled" SMS operator (atomic units): 

Q=, — - — r +Ap x -p 2 =J2Qk, (A.l) 
|ri-r 2 | ^ 

where A is the scaling factor, p = — zV is electron momentum, and 

Qk = 2k + lrk l 1 Y k (n 1 )Y k (n 2 ) + A Pl ■ p 2 5 kl . (A.2) 
We use the following form for the single-electron wave function 

I (A.3) 

Here a = 1/137.036 is the fine structure constant, and Cl(n)ji m = —{a ■ n)fi(n) J 7 m . 
The matrix element of operator (jA.2|) with wave functions (|A.3|) has the form 

{^x{v x )^ 2 {v 2 )\Q k \^{v x )^{v 2 )) = C k (R k - APi4i), (A.4) 

where the angular factor C k is the same for both operators 

C = (-1)3+^1+^2 ( Jl k 32 k 34 ^ 

k ^ ' I —TTli 1 m 3 J V —177-2 —Q m 4 ) 



x(-l)^+^»+w+Y(2j 1 + l)(2j 2 + l)(2j 3 + l)(2j 4 + 1) (A.5) 

X (f -I o)(f -\ o)^ h + h + k ^ + h + k ^ 

I 1, if x is even, 



ax 

R k is radial Coulomb integral 



0, if x is odd 



POO /r>k 

Rk = / ^r(/i(ri)/ 3 (ri) + a 2 g 1 (r 1 )g 3 (r 1 ))(f 2 (r 2 )f 4 (r 2 ) + a 2 g 2 (r 2 )g 4 (r 2 ))d ri dr 2 , (A.6) 
Jo r> 

while Pi is radial matrix element of the SMS operator 

Pi = PnP24, 

Pab = A ab 5 Ub+ i + B ab 5 Ub -i, (A. 7) 

A-ab = Jo fa(jp — l -f)fbdr, 

B ab = J oo f a (i + ^)fbdr. 
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